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SUMMARY 


A necessary parameter to be investigated in the analysis of space 
missions is the total assembled mass in Earth orbit, as it is a strong indi- 
cator in comparisons of the performance of various propulsion systems. The 
computer code described in this paper was designed to compute the mass in 
Earth orbit for round-trip stopover missions, flyby missions, and orbiter 
missions with a single-stage high-thrust system for each propulsive phase. 
The actual gravity losses suffered during each propulsive phase are computed 
by complete integration of the equations of motion, thereby permitting the 
propulsion system to be optimized. The program allows the input of planet 
ephemeris data so that it may be used for any planet. The following options 
are available to the user. 

The number of engines may be specified or optimized 
for each propulsive phase of the mission. 

Planet capture may be performed by aerobraking 
or propulsive maneuvers. 

A specific number of engines/tanks may be 
specified. 

The engine burn time may be constrained to a 
desired value. 

Midcourse correction penalty may be specified. 

Time dependent life support may be included for 
manned missions. 

The inert weight fractions may be specified as 
constants or computed by the program. 

An engine clustering penalty may be specified. 

Arrival or departure or both from elliptical orbits 
may be specified. 

Thrust may be initiated or terminated at points other 
than per laps is. 



Perturbing of the parking orbit may be investigated, 
i.e., by using an apoapsis kick. 

Each phase of a mission may be run independently. 

Thrust vector orientation other than tangential may 
be specified. 

Examples provided include the complete input and output for a round-trip 
mission, a flyby mission, and an orbiter mission. 


INTRODUCTION 


The analytical computation of the assembled mass in Earth orbit for space 
missions is time consuming and difficult. Approximations are used in such 
computations for the gravity losses suffered during propulsive phases, the 
inert weight scaling laws, etc. Therefore, a computer code employing numeri- 
cal integration of the equations of motion, optimization procedures on several 
variables, and iterative procedures for inert calculations has been developed 
that enables the user to compute accurately the mass in Earth orbit for round- 
trip missions, orbiter missions, and flyby missions. 

The flyby type mission consists of one single propulsive phase, Earth 
departure, and is assumed to place a specified weight on a heliocentric flyby 
trajectory. The orbiter type mission consists of two phases, Earth departure 
and planet capture. Both phases may be propulsive or the capture phase may 
be performed by an aerobraking maneuver. A midcourse correction fuel penalty 
may be provided. 

The round-trip stopover mission has three phases, Earth departure, planet 
capture, and planet departure, and includes a provision for an excursion mod- 
ule for landing on the target planet. The excursion module weight must 
include the propulsion necessary to descend to and ascend from the planet and 
any staged weight left at the planet (i.e., small probes). The Earth return 
weight would include the mission module Earth entry vehicle and time dependent 
expendables (e.g., life support). Any propulsion fuel consumed in midcourse 
can be included. This computer code has been written in Fortran IV and has 
been used on a 709^ computer. The fact that Fortran IV was used as a program 
language should make it applicable to other types of computing equipment . 

Required program input, output, example problems and the program listings 
are included here as appendixes A-D, respectively. 
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ANALYSIS 


Method of Solution 


The procedures for computing total assembled mass in Earth orbit for 
round trip stopover miss ions , orbiter missions , and flyby missions are out- 
lined in this paper. The program performs the calculations for the round trip 
missions in three distinct phases (i.e., phase I, planet departure, phase 2, 
planet capture, and phase 3> Earth departure), starting with the Earth return 
weight and working in reverse order through the total missions. The orbiter 
mission consists only of phases 2 and 3 and a flyby mission only of phase 3 . 
The inputs that the user must provide include propulsion system characteris- 
tics (e.g., engine weight, engine thrust), trajectory data (e.g., Voo, orbit 
altitudes), mission module weights, staged weights at planet and time depen- 
dent expendables (e.g., life support). A complete listing of all required and 
optional input is included in appendix A. 

The actual gravity losses suffered during each propulsive phase of a 
mission are computed by numerical integration of the equation of motion in 
order to determine the optimum ratio of thrust to weight of each stage, the 
inert fractions, propellant consumption, and the engine burn times. The pro- 
grams also have an option that allows the mass in Earth orbit to be computed 
by the closed form ballistic equations rather than integration of the equa- 
tions of motion, and includes an approximation for the gravity losses in a 
separate subroutine. Other available options are described in the section 
Program Options . 

This computer code is capable of computing four separate types of 
missions. The mission type must be specified by the input parameter MODE as 
follows : 

MODE = 1 Round trip stopover missions 
MODE = 2 Orbiter missions 

MODE = 3 Flyby missions 

MODE = 4 Only the planet capture phase is computed 

The method of solution described below is for a round trip stopover mission 
with a specified Earth return module which contains an Earth entry vehicle. 
Starting with the Earth return module the program works through the mission 
profile in reverse order in the following manner. The planet departure pay- 
load is set equal to the weight of the Earth return module plus the midcourse 
fuel and time dependent expendables . An initial estimate of the ratio of 
total engine thrust to total vehicle weight is computed by the program for the 
planet departure phase. This initial estimate is computed from the closed- 
form ballistic equation with no gravity losses. The trajectory is then inte- 
grated using the estimated t/Wi until the specified hyperbolic excess veloc- 
ity for planet departure is obtained, at which time the propulsive factor P 
is computed . 
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p 


1.0 + (we/wl) 

MFST (l + A) - A 


where 

WE engine weight 

WL Earth return weight + time dependent expendables + midcourse fuel 

A inert fraction (either input or allowed to be computed) 

MFST propellant consumed during powered phase 
and 

MFST = 1.0 - T /^-* - 
Isp 

t/w engine thrust/gross weight of vehicle 

t burn time* sec 

Isp specific impulse, sec 

This propulsion factor is then multiplied by the returned payload weight to 
arrive at a total gross weight of the stage based on the previous integration. 
Forming the t/w 2 based on the new found gross weight of the stage, the pro- 
gram compares the initial estimate for t/Wi to the t/w 2 computed to see 
whether the two values agree within a desired tolerance. If the values dis- 
agree by more than the desired tolerance, the trajectory is rerun with t/w 2 
and the iteration continues until convergence is obtained. The value of the 
gross weight, WG = P * WL, is then taken as the gross weight of the planet 
departure stage. The ratio of propellant weight to gross weight WP/WG is 
given by WP/WG = 1.0 - MFST and the weight of propellant is computed from 
WP = WP/WG * WG. The weight of the inerts (less the engines) is found from 
WT = WP * A where A is the inert fraction computed from the scaling law 
equations or is an input constant. 

The payload weight of the capture phase of the mission is set equal to 
the gross weight of the planet departure phase plus the weight of any payload 
left at the planet. 


WIA = WG + WEM 

where WEM is the weight of staged payload (e.g., planet excursion module). 
An initial weight of propellant is assumed by the program and the initial 
estimate for the WTG, inert weight, is made, using the inert fraction com- 
puted from the scaling laws. Using these initial estimates, an estimate of 
the total gross weight of this stage less the propellant is given by 


k 


WGL = WL1 + WE + WTG 


where 


WE weight of engines 

WTG initial estimate of inert weight (less engine) 

WL1 gross weight of departure stage + staged payload. 

The value of WG2 is used to form the thrust-to-weight ratio to start the 
integration procedure. The equations of motion are integrated using negative 
time until the desired hyperbolic excess velocity, is obtained for the 
capture phase. The use of -t in the integration allows the equation 

MFST = 1.0 - (T/W * t/lsp) 

to add propellant to the vehicle as the trajectory is computed in reverse 
order starting from the low capture orbit. When the desired is obtained 

the propulsive factor P is computed 


1 + (we/wll) 

(l.o/MFST) (1.0 + A) - A 


The product P * WLl then gives the total weight, WBO, of the capture stage 
prior to the start of the capture maneuver. The total weight of propellant 
used during the capture phase is 


WFC = WBO - WG2 

and the weight of the inerts (less the engine), WTG, is 

WTC = WPC * A 

where A is the inert fraction. This value of WTG is compared to the ini- 
tial guess WTG to see whether the values are within the desired tolerance. 
If the tolerance is met, the gross weight of the capture phase of the mission 
is the WBO found. If the tolerance is not met, WTC is used as the next 
guess of WTG and the iteration continues until convergence is obtained. 

When the third phase of the mission, Earth departure, is initiated, the 
gross weight of the capture phase of the mission is used as a base weight and 
is increased as follows . The payload of the third phase is 

WLE = WBO + WM •+■ WLS 


where 


5 



WBO gross weight of capture phase 

WM weight of midcourse fuel 

WLS weight of time dependent expendables (e.g., life support) 

This value of WLE is then used to form the thrust /weight needed to start the 
integration of the equation of motion for the Earth departure phase. The 
iteration procedure described in the planet departure phase of the mission is 
the same employed for the Earth departure and all values of inert weights, 
propellant weights and gross weights are found by the same method. 


Program Options 

There are several options in this computer code that are available to the 
user. They are listed here and explained further in the section devoted to 
the computer input. 

(1) A limit on engine operating time has been set in the program to 1800 
sec. This may be changed by reading in TLIMIT equal to the desired operating 
time limit in seconds. When this limit is exceeded, the number of engines 
used is increased by one until the operating time constraint is met. 

(2) The number of engines used per propulsive phase may be fixed or 
optimized at the discretion of the user. When a fixed number of engines is to 
be used (no optimization), set SEARCH = 0 and N^El, N$E2, N^E3 equal to 
desired number of engines. (Here and in the following 1, 2, and 3 refer to 
planet departure, planet capture; and Earth departure maneuvers, respectively.) 
The computer code will not change the number of engines specified unless the 
operating time limit is exceeded; therefore, set TLIMIT to a large positive 
number if no operating limit is desired. The user may set SEARCH = 1 and the 
engine optimization procedure is followed; that is, the number of engines 
specified by N$E1, N0E2, N$E3 is used as a starting base and the number of 
engines is reduced by one until a minimum gross weight is obtained. The opti- 
mum number of engines is then output with the associated gross weight. The 
operating time limit will also be a factor and should be set to the desired 
limit when optimizing the number of engines. 

( 3 ) Either a specific engine thrust level may be used as noted in (2) 

above, or the engine may be " rubberized 3 * * * * * * * 11 as follows: Estimate a total thrust 

level needed for the phase of the mission being computed. Divide this thrust 

level by the value of N^El, H0E2, or N</E3 pertinent to the mission phase and 

set ET1, ET2, or ET3 equal to the thrust thus obtained. Also, set SEARCH = 1, 

allowing the program to optimize the number of engines . If the output lists 

the number of engines as 1, choose a lower thrust per engine and rerun the 
problem. If the number of engines is equal to the input value, choose a 
higher thrust per engine and repeat the run. When the output number of 

engines is greater than 1 and less than the input value, the number of 

"rubberized” engines has been optimized. 
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(4) Aerobraking for the capture phase of a mission may he used by setting 
IAER0 = 2 and assigning a heat shield penalty to the incoming vehicle weight. 
If a penalty of 15 percent of the vehicle weight was desired, the user would 
input AER0P = 1.15- 

(5) The inert fraction (less the engine) may be an input or may be 
computed in the program from the following equations which were developed by 
United Aircraft Corp. under contract NAS 2-2928. Tankage structure weight 
WS : 


WS 


A*TH° ' 1 0-9 

0 O.533 


+ Ki TN° ' 1 


where 

A constant 

a specific gravity of propellant 

Ki constant for fixed weights 

WP propellant weight 

TN number of tanks 


The number of tanks is set initially at 1.0. If the user specifies the number 
of engines per tank, NET1, then the program will determine the optimum number 
of tanks as well as the optimum number of engines. If a value for NET1, NET2, 
or NET3 is not read as input, the program will compute the mission using one 
single tank for each phase that has not specified a particular number of 
engines per tank. 


(6) Boil off and thermoinsulation for cryogenic tanks 


WBO + Wp 



where 

B constant in kms system 

WP weight of usable propellant 

a specific gravity of propellants 

t time of exposure, days 

L latent heat of vaporization, kms system 

AT temperature difference across insulation thickness 
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(7) To avoid integration of micrometeoroid impacts with respect to 
exposure time in estimating weight of meteoroid protection, and at the same 
time include the effects of the relatively large fluxes of debris in the 
asteroid belt, the following approximate equation has been used: 


WM = 


C 



1/4 


The choice of the values assigned to the constant, C, should distinguish 
between missions to regions of the solar system inside and within or beyond 
the asteroid belt. The following values of C are suggested. 

Parameter C 

0 to 1.8 A.U. 1.8 A.U. & beyond 


High 0.06 0.10 

Nominal .03 *05 

Low .015 .025 

These values of C are for single-sheet thicknesses, with time t measured 
in hours. Inasmuch as the tank structure WS affords some protection against 
meteoroid penetration, only an additional weight, AWM, need to be added to 
the tank structure. If the calculated meteoroid protection, WM, is less than 
the tank structure, WS, no additional weight is added to the tank. If WM is 
greater than WS but less than 1.33 WS, additional weight, AWM = WM - WS, is 
added to the single tank sheet. If the required meteoroid protection is 
greater than 1.33 WS, a second tank sheet is used to construct a Whipple 
bumper with the additional weight AWM = WS/3. 

(8) For arrival at or departure from elliptical orbits the periapsis and 
apoapsis radii (KPP and RAP, respectively) are specified. If the true anomaly 
is not specified, the program will initiate thrust at the periapsis. To ter- 
minate or initiate thrust at points along the orbit other than periapsis TLT1, 
TLT2, and TLT3 must be specified equal to the desired true anomaly. 

(9) The user may also perturb the parking orbits by such maneuvers as an 
apogee kick by placing KICK1, KICK2, or KICK3 equal to the change in periapsis 
radius desired at the respective orbit. 

(10) If fixed inert fractions rather than calculated inerts are desired, 
simply place AC^MP = 1 and input the desired inert fractions Al, A2, and A3. 

(11) A fixed Earth departure stage, such as a Saturn SII, may be treated 
by placing AT3 =0.0 and WK3 = 37500.0 which is approximately the fixed tank 
inerts on an SII stage. With AC0MP = 2, this same procedure may be used at 
any phase of the mission. 

(12) A useful option at the user's command is the ability to operate the 
program at two different levels of accuracy. With LEVEL = 1, all propulsive 
phases of the mission are computed in closed form with an estimate of the 
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gravity losses, which is an empirical curve fit to the results of reference 1 . 
With LEVEL = 2, all propulsive phases are computed with a complete integration 
of the equations of motions. The first level is useful for scanning many 
missions because of the minimal machine time. The second level, although 
somewhat longer in computer time, allows a more exacting analysis of a mission. 
LEVEL = 2 is built into the code and will perform the integrations if it is 
not changed by the input LEVEL = 1 as noted above. 

(13) A clustering penalty of 10 percent per additional engine has been 
included in the program but may be changed by input. If for example a 12- 
percent penalty were to be used, set PEN = 0.12 in the mission data input. 


General Comments 


Stacked cases .- Cases may be stacked by simply placing an asterisk card 
between each set of input. Only those input quantities which are to be 
changed need be added after the asterisk. If no input is read for the inert 
calculation, two ( 2 ) asterisk cards must follow the stacked mission data. 

Perturbed inputs .- The program will utilize the last input card 
encountered for a given parameter in its analysis. Hence, if a case requires 
the changing of one or more input cards from a prior case, one can simply 
place the desired changes on a card at the end of the pertinent input and it 
will utilize these latter values in the analysis. Thus, whole problem decks 
can be duplicated for stacked cases and perturbations made using supplementary 
cards without having to repunch the original problem deck. 

All the variables in the Input Parameter List, appendix A, have been 
given mode numbers in the left-hand column, i.e., 1, 2, 3 , k. These refer to 
the type mission requiring that particular variable. The four modes are as 
follows : 


M0DE = 1 
M0DE = 2 
M^DE = 3 
M0DE = 4 


Round-trip mission (all variables marked l) 
Orbiter mission (all variables marked 2) 

Flyby mission (all variables marked 3) 

The capture phase only is computed (all variables 
marked 4 ) . 


It should be noted that some of the variables are needed for all four mission 
modes while some are only necessary for one particular mission mode. Mode 4 
may be used for investigation of propulsion system characteristics for the 
capture phase only. If a survey of only the departure phase is desired, the 
user may set MODE = 3> and GSE, RGE, RPE and RAE equal to the departure 
planet's respective values. 


National Aeronautics and Space Administration 

Moffett Field, Calif., 94035, Nov. 7, 1967 

789 - 4 O-OI-O 8 -OO-I 5 
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APPENDIX A 


INPUT PARAMETERS 


Mode Variable name 


Description 


1 




EWl 

1, 

2, 

4 


EW2 

1, 

2, 

3 


EW3 

1 




ET1 

1, 

2, 

4 


ET2 

1, 

2, 

3 


ET3 

1 




N0E1 

1, 

2, 

4 


N0E2 

1, 

2, 

3 


N0E3 

1 




11 

1, 

2, 

4 


12 

1, 

2, 

3 


13 

1 




A1 

1, 

2, 

4 


A2 

1, 

2, 

3 


A3 

1 




V8C01 

1, 

2, 

4 


V8C02 

1, 

2, 

3 


V8C03 

1 




MCPF1 

1, 

2, 

3 


MCPF 

1, 

2, 

3, 

4 

WPL 

1 




WEM 

1, 

2, 

3, 

4 

GSE 

1, 

2, 

3> 

4 

RAE 

1, 

2, 

3, 

4 

RPE 

1, 

2, 

3> 

4 

RGE 

1, 

2, 

4 


RAP 

1, 

2, 

4 


RPP 

1; 

2, 

4 


RGP 

1, 

2, 

4 


GSP 

1 




FELL0W 

1 




LEGT1 

1 




LEGT2 

1 




LSD 

1, 

2, 

4 


AER0P 

1, 

2, 

3, 

4 

SEARCH 


Weight /engine, lb, planet depart 
Weight/engine, lb, planet capture 
Weight /engine, lb, Earth depart 
Th rust/ eng ine , lb, planet depart 
Thrust/ engine, lb, planet capture 
Thrust/ eng ine, lb, Earth depart 
Number of engines , planet depart 
Number of engines , planet capture 
Number of engines , Earth depart 
Isp, sec, planet depart 
Isp, sec, planet capture 
Isp, sec, Earth depart 

Inert fraction (less engine), planet depart 
Inert fraction (less engine), planet capture 
Inert fraction (less engine), Earth depart 
Voq, km/sec, planet depart 
V , km/ sec , planet capture 
Voo, km/sec. Earth depart 
l.XX - midcourse correction penalty, 
inbound leg 

l.XX - midcourse correction penalty, 
outbound leg 

Weight of MM + EEM, etc., lb 
Weight of MEM, lb 
"g” at Earth, km/sec 2 

Apogee of departure orbit at Earth, © radii 
Perigee of departure orbit at Earth, © radii 
Radius of Earth, km 

Apogee of capture orbit at planet, planet 
radii 

Perigee of capture orbit at planet, planet 
radii 

Radius of destination planet, km 
"g" at destination planet, km/sec 2 
Number of men 
Inbound leg time , days 
Outbound leg time , days 
Life support requirements, lb/man-day 
l.XX - penalty for aerobraking at destina- 
tion planet 

0 - will not optimize number of engines but 
use input values 

1 - will optimize number of engines 
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1, 2, 3> 4 AC0MP 1 - will use fixed inerts 

2 - will calculate inerts using inert 
subroutine 

1, 2, 4 IA.ER0 1 - propulsive capture at destination 

planet 

2 - aerobraking at destination planet 
1, 2, 3, 4 PEN Engine clustering penalty 

1, 2, 3, 4 TLIMIT Maximum allowable engine bum time, sec 

1 , 2, 3, 4 RET1 Number of engines/tank, planet departure, XX. X 

1, 2, 4 NET2 Number of engines/tank, planet capture, XX. X 

1 , 2, 3 NET3 Number of engines/tank. Earth departure, XX. X 

1, 2, 3, 4 M0DE The type mission to be computed 

1, 2, 3, 4 LE VEL 1 - propulsive phases computed in closed 

form with approximation for gravity loss 
2 - propulsive phases computed by numerical 
integration of equation of motion 

1 KICK1 Desired % change in periapsis radius, planet 

depart 

1, 2, 4 KICK2 Desired % change in periapsis radius, planet 

capture 

1, 2, 3 KICK3 Desired % change in periapsis radius, Earth 

depart 

1 TLT1 True anomaly at initiation of thrust, planet 

depart, deg 

1, 2, 4 TLT2 True anomaly at initiation of thrust, planet 

capture, deg 

1, 2, 3 TLT3 True anomaly at initiation of thrust, Earth 

depart, deg 

1, 2, 3, 4 * An asterisk in any column (1-72) will termi- 

nate the reading of the above mission data. 

The following input quantities are used in the inert calculations . This 
data must follow the above mission data. If AC^MP = 1 in input above, the 
constant inert fractions Al, A2, A3 above are used, and the following input 
may be deleted. 

1 SIG1 Propellant specific gravity, planet depart 

1, 2, 4 SIG2 Propellant specific gravity, planet capture 

1, 2, 3 SIG3 Propellant specific gravity. Earth depart 

1 ATI Constant (A) in tankage equation, planet 

depart 

1, 2, 4 AT2 Constant (A) in tankage equation, planet 

capture 

1, 2, 3 AT3 Constant (A) in tankage equation, Earth 

depart 

1 B01 Constant (b) in boil off equation, planet 

depart 

1, 2, 4 B02 Constant (B) in boil off equation, planet 

capture 

1, 2, 3 B03 Constant (B) in boil off equation. Earth 

depart 
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1 


CM1 

1, 2, 4 CM2 

1, 2, 3 CM3 

1 TEMPI 

1, 2, 4 TEMP2 

1, 2, 3 TEMP3 

1 HEAT1 


2 , 4 
1 , 2 , 3 
1 

1 , 2 , 4 
1 , 2 , 3 
1 

1 , 2 , 4 

1 , 2., 3 
1 , 2 , 3 , 4 


HEAT2 

HEAT3 

TEX1 

TEX2 

TEX3 

WK1 

WK2 

WK3 

* 


Constant (c) in meteoroid equation, planet 
depart 

Constant (c) in meteoroid equation, planet 
capture 

Constant (C) in meteoroid equation, Earth 
depart 

AT in boil -off equation, °K, planet depart 
AT in boil-off equation, °K, planet capture 
AT in boil-off equation, °K, Earth depart 
Propellant latent heat of vaporization (l) 
in boil-off equation, kcal/kg, planet 
depart 

Propellant latent heat of vaporization (L) in 
boil-off equation, kcal/kg, planet capture 
Propellant latent heat of vaporization (L) in 
boil-off equation, kcal/kg. Earth depart 
Exposure time (t) boil-off and meteoroid 
equation, days, planet depart 
Exposure time (t) boil-off and meteoroid 
equation, days, planet capture 
Exposure time (t) boil-off and meteoroid 
equation, days. Earth depart 
Constant K in tankage equation, kg, 
planet depart 

Constant K in tankage equation, kg, planet 
capture 

Constant K in tankage equation, kg, Earth 
depart 

An asterisk in any column (l-72) will terni- 
nate the reading of the inert data and 
execution will begin 
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APPENDIX B 


OUTPUT PARAMETERS 


The following output will appear for each phase of the mission computed. 


Variable name 

Description 

Engines 

Optimum number of engines 

Thrust 

Total thrust , lb 

Engine weight 

Total engine weight, lb 

VINF 

Vo,, km/sec 

VGL 

Gravity losses, km/sec 

DEL V 

AV including gravity losses, 
km/ sec 

Burn time 

Stage bum time, sec 

Isp 

Stage Isp, sec 

A 

Inert fraction (exclusive of 
engine wt) = Wj/Wp 

MUL 

Payload fraction = Wp/WQ 

Tanks 

Number of tanks for that stage 

Payload 

Total stage payload, lb 

W TANK 

Weight of all inerts except 
engine (W«j + W^q + Wj + Wjyjjfl), 
lb 

W FUEL 

Weight of propellant used, lb 

GR0SS WT 

Gross weight, lb 

P 

Propulsion factor = l/nL 

t/w 

Stage initial thrust to weight 

MID Fuel 

Midcourse correction penalty, lb 

LIFE WT 

Life support weight , lb 

THETA 

Angle between original line of 
apsides and departure 
asymptote 

EX.M0D 

Excursion module weight, lb 

Tank 

Wp/Wp 

Boil off and insulation 

(Wpo + Wj)/Wp 

Meteoroid shield 

W MM' W P 
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APPENDIX C 


EXAMPLE PROBLEMS 


Included in this section are sample runs of a manned roundtrip, an 
unmanned orb iter, an unmanned flyby, and a capture phase only. It is noted 
that an asterisk card should appear between the mission input data and the 
system inert parameter data. 

In the first example, nuclear propulsion was used at all phases; however, 
the thrust level of an individual engine was different. Planet departure and 
capture was performed with a small 75; 000 lb thrust engine while the Earth 
departure phase used a 200,000 lb thrust engine. The input SEARCH = 1 was 
used allowing the number of engines to be optimized. The operating time limit 
is set to 1,000 sec, and the number of engines per tank NET1 = NET2 = NET3 = 
2 . 0 . 


The program increased the number of engines to two for the planet 
departure phase and to three for the planet capture phase because fewer 
engines exceeded the operating time constraint. 

The optimization indicates three engines, for the Earth departure phase 
produced the minimum total gross weight while remaining within the operating 
time limit. The separate inert fraction for tankage, boil off and insulation, 
and meteoroid shielding are output for the user 1 s reference . Note that the 
boil off and insulation fraction are zero for the Earth departure phase 
because TEX3 =0.0 (the exposure time) was assumed. 

Example II is for a Mercury orbiter mission. Note that the planet 
parameters RAP, RPP, RGP, and GSP have been changed to those for Mercury. The 
parameter from the input list with a mode 2 designation only have been sup- 
plied. The engine operating time limit was not specified and the 1800 sec 
built in was used. This example will have one tank only since the number of 
engines per tank NET2, NET3 was not read in. The output consists only of two 
phases printed because no planet departure phase was computed. The input 
value M0DE = 2 was required to accomplish this orbiter mission. 

Example III is a flyby to Mercury. The value M0DE = 3 is read in and 
only those input parameters for mode 3 are needed. The output as noted will 
consist of only an Earth departure phase. The inert fractions printed show a 
value for the tankage while the boil off and insulation and meteoroid shield- 
ing are zero. The boil off and insulation are zero as TEX3 (exposure time) 
was set to 0.0 and the meteoroid shielding is zero, indicating the tankage 
will afford the needed meteoroid shielding based on the scaling law. TEX3 
may be input as the assembly time in Earth orbit for manned mission or the 
parking orbit time for unmanned missions. 
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EXAMPLE OF INPUT FOR ROUND-TRIP MISSION (MODE=l) 
VT I T =60H SAMPLE MARS STOPOVER MISSION , 

MODE= 1 t 

E W 1= 17000.0, E T 1=75000 « 0 , N0E1=1, 

E W2= 17000 . 0 , ET2=75000.0, N0E2=1* 

EW3= 28700.0, E T 3= 200000 .0 , NOE 3 = 6, 

I 1 = 850.0, I 2=8 50.0, 13=850.0, 

A 1= . 25 , A2= . 2 5 » A3=.25, 

V SCO 1 = 6.92, V8C02=3. 59, V8C03=2.86 5, 

MCPF = 1.02, MCPF 1=1.02, 

W PL =92000.0, HEM= 134000.0, 

RAE = 1 .066 » RPE= 1.066, 

R AP = 1 . 265 , RGP = 3410.0, RPP=1.265, GSP=. 00369, 

SF ARC H= 1 , ACHMP = 2 , 

FFLL0W=8.0, LEGT 1=260.0, LFGT2=180.0, LSD=3.75, 

PEM=. 1, 

NET 1 = 2.0, NFT2 = 2.0, NET3 = 2.0, 

T L I M I T = 1000 . 0 , 

INPUT FOR INFRT ROUTINE 
S IG 1= . 077 , S IG2= • 077 , SIG3=.077, 

AT 1= . 10 , AT 2= • 10 , AT 3= . 10 , 

B01= .034, B02=.034, B03=.034, 

CM 1= . 06 , CM 2= .06 , CM3=.06, 

TEMP1=156.0, TE MP 2= 1 56 . 0 , TFMP3=156.0, 

HE AT 1= 108 . 1 , HE AT2= 108 . 1 , HEAT3=108.1, 

TEX1=210.0, TEX2=180.0, TEX3=0.0, 

WK 1 = 500.0, WK2 = 500 .0 , WK3=500.0, 


SAMPLE MARS STOPOVER MISSION 

DEPARTURE FROM PLANET 

ENGINES THRUST ENG. WT VINE VGL DEL.V BURN TIME ISP A MUL TANKS 

2 150000.0 37400.0 6.920 0.0712 5.1502 855.75 850.0 0.262 0.3059 1.0 

PAY LOAD W TANK W FUEL GROSS WT P T/W MID FUEL LIFE WT THETAP 

101640.0 40061.2 153119.0 332220.2 3.269 0.4515 1840.0 7800.0 114.6 

INERT FRACTION A CONSISTS OF TA NK = 0 . 1 36 192 BOIL OFF + I NSUL AT 1 0N = 0 . 080044 METEROID SH I ELD = 0. 04 540 


CAPTURE PHASE 


ENGINES 

3 

THRUST 

225000.0 

ENG. WT 
56100.0 

VINF VGL 

3.590 0.0265 

DEL.V BURN TIME ISP A 

2.5982 796.66 850.0 0.258 

MUL TANKS 

0.5918 2.0 

PAY LOAD 
466220.2 

W TANK 
54540.2 

W FUEL 
211172.9 

GROSS WT 
787763.9 

P T/W EX MOD THETAP 

1.690 0.2856 134000.0 -257.9 



INERT FRACTION A CONSISTS OF T A NK = 0 . 14401 8 BOIL OFF + I NSUL AT I0N = 0 .066249 METEROID SH I EL 0 = 0. 04 8 0 1 

DEPARTURE FROM EARTH 

ENGINES THRUST ENG. WT V1NF VGL DEL.V BURN TIME ISP A MUL TANKS 

3 600000.0 94710.0 2.865 0.0558 3.6006 736.33 850.0 0.126 0.5415 2.0 

PAY LOAD W TANK W FUEL GROSS WT P T/W MID FUEL LIFE WT THETAP 

808919.2 66260.7 523990.0 1493879.9 1.847 0.4016 15755.3 5400.0 174.4 

INERT FRACTION A CONSISTS OF T A NK=0 . 1 26454 BOIL OFF + I NSUL AT I0N=0 . 000000 MFTEROIP SHI ELD=0. OUOOO 
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EXAMPLE OF INPUT FOR ORBITER MISSION (MODE=2> 
VT I T =6 OH SAMPLE MERCURY ORBITER MISSION , 

M0DE=2, 

EW2=312.0, ET2=25000.0, N0E2=6, 

EW3=28700.0, E T3=200000 . 0, N0E3=1, 

12=450 .0 , 1 3=850 *0 * 

V SCO 2 = 8 . 04 f V 8C03= 1 1 .47 , 

MCPF=1.02, 

WPL= 10000 .0 t 
RAE=1.036, RPE=1.036, 

RAP= 1.041, RPP= 1,041 , RGP=2420.0 » GSP*. 00372, 

SE ARCH= 1 , AC0MP=2, 

• 

INPUT FOR INERT SUBROUTINE 
S IG2= .077 , S IG3= ,077 , 

AT 2= . 10 , AT3= . 10 , 

B02= .034, B03=,034, 

CM2= ,06 , CM3=.06, 

TEMP2= 156,0 * TEMP 3=156,0, 

HE AT2= 108 .It HEAT3= 108 . 1 , 

TEX 2= 130.0* TEX3=0 ,0 , 

WK 2 = 5 00 . 0 , WK3= 500 « 0 , 

S IG2=, 33, 

TEMP2= 104 .0 * 

HEAT2=60,5* 


SAMPLE MERCURY ORBITER MISSION 

CAPTURE PHASE 

ENGINES THRUST ENG. WT VINF VGL OEL.V BURN TIME ISP A MUL TANKS 

2 50000.0 686.4 8.040 0.0614 6.1726 487.98 450.0 0.131 0.1385 1.0 

PAY LOAD W TANK W FUEL GROSS WT P T/W EX MOO THETAP 

10000.0 7141.7 54404.0 72180.8 7.218 0.6927 0.0 -271.3 

INERT FRACTION A CONSISTS OF TANK=0. 086059 BOIL OFF + INSULAT ION=0. 036615 MET EROI D SHI ELD=0. 00860 

DEPARTURE FROM EARTH 

ENGINES THRUST ENG. WT VINF VGL DEL . V BURN TIME ISP A MUL TANKS 

1 200000.0 28700.0 11.470 0.3081 8.4223 981.14 850.0 0.128 0.2036 1.0 

PAY LOAD W TANK W FUEL GROSS WT P T/W MID FUEL LIFE WT THETAP 

74010.6 29652.0 231194.2 363556.8 4.912 0.5501 1451.2 0.0 138.5 

INERT FRACTION A CONSISTS OF TANK=0. 1282 56 BOIL OFF + INSULAT I ON =0.000000 MET EROI D SHI ELD=0. 00000 
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EXAMPLE OF INPUT FOR FLY-BY MISSION ( M0DE=3 ) 
VT IT =60H SAMPLE MERCURY FLYBY MISSION , 

MODE =3 , 

EW3= 17000.0, ET3=75000.0, NOE 3=1, 

13=850.0, 

V0CO3= 1 1 .47 , 

MCPF= I . 02 , 

WPL= 10000 .0 , 

R AE = 1. 036 , RPE = 1 • 036 , 

RAP-1.041, RPP= 1.041 , RGP = 2420 .0 , GSP=. 00372, 

S E ARCH= 1 , ACOMP= 2 , 

* 

INPUT FOR INERT SUBROUTINE (P002IN) 

S I G3= .077 , 

AT3=. 10, 

B03 = .034 , 

C M3= . 06 , 

TEMP3=156.0, 

HE AT3= 108 . 1 , 

TE X3 = 0 .0 , 

WK3=500 .0 , 


SAMPLE MERCURY FLYBY MISSION 


DEPARTURE FROM EARTH 


ENGINES 

1 


THRUST 

75000.0 


ENG. WT 
17000.0 


VINF VGL DEL . V 

11.470 0.1839 8.2981 


BURN TIME 
720.58 


ISP A MUL TANKS 

850.0 0.158 0.1013 1 


PAY LOAD W TANK 

10200.0 10020.7 


W FUEL GROSS WT 
63500.7 100721.3 


P T/W MID FUEL 

9.875 0.7446 200 .0 


LIFE WT THETAP 

0.0 131.2 


INERT FRACTION A CONSISTS OF TANK = 0. 157804 BOIL OFF + INSUL AT 1 0N=0 .000000 METEROID SHI ELD=0. 00000 


EXAMPLE OF INPUT FOR CAPTURE PHASE <M0DE*4) 
VT IT =60H EXAMPLE CAPTURE PHASE ONLY 
MODENA, 

EW2 = 3 12 .0 , ET2=25000.0, N0E2=6, 

12=450.0, 

V8C02=8.04, 

W PL= 10000*0 , 

RAE = 1 . 0 36 , RPE = 1.0 36, 

R AP= 1.041, RPP=1.041, RGP=2420.0, GSP=. 00372* 

SE ARCH= 1 , AC0MP*2, 

* 

INPUT FOR INERT SUBROUTINE (P007IN) 

AT2=. 10, 

B02=.034, 

CM2=.06, 

TE X2= 1 30 ,0 , 

WK2=500.0 , 

S I G2* • 33 * 

T EMP2= 104.0, 

HE AT2=60 .5 , 


EXAMPLE CAPTURE PHASE ONLY 


CAPTURE PHASE 


ENGINES 

THRUST 

ENG. WT 

VINE 

VGL 

DEL .V 

BURN TIME 

ISP 

A 

MUL TANKS 

2 

50000.0 

686.4 

8.040 

0.0614 

6.1726 

487.98 

450.0 

0.131 

0.1385 


PAY LOAD W TANK W FUEL GROSS WT P T/W EX MOO THETAP 

10000.0 7141.7 54404.0 72180.8 7.218 0.6927 0.0 -271.3 


INERT FRACTION A CONSISTS OF 


TANK=0. 086059 


BOIL OFF + INSULATION=0. 036615 


METEROI D SHIELD=0. 00860 



APPENDIX D - PROGRAM LISTINGS 


MAIN CONTROL PROGRAM FOR REQUIRED MASS IN ORBIT COMPUTATION I 

COMMON ALTO,PSI 1,1 12,GE,GM,ALT,T(70) ,PSI2,K,GS 2 

COMMON/ BO XGL/ DVI , F W 12, V SCO , VSU , DVGL ,L EVEL , T ST , RG , R LTS , AS 3 

COMMON/ BOX 1 /A 1, A2, A3, A 12, I ABC , ACOM P , $ I G 1 , S I G2 , S I G 3 , ATI , BO 1 , CM1 , AT2 A 

1 , 602 , CM2 , AT3 , BOB , CM3 , TE M P 1 , T EM P 2 , TEM P3 , HEAT 1 » HE AT2 , HE AT3 , T EX 1 , TEX2 5 

2, TEX3,WK1,WK2, WK3,TN 6 

DIMENSION V T I T ( 18)»VIN8(20) ,MF8( 10) ,TH018( 10)»TVIN8(10),TMF8(10)»T 7 

1TH0L8I 10),T99( 10),TT99( 10) 8 

REAL 111,1 12,MF12,ML10,MLll,ML12tMLl,K,KICK ,MUL ,MFST,MF8,I 1,12,13, 9 

1MCPF,MCFW,KICK1,KICK2,KICK3,MCPF1,MCFW1,LEGT1,LEGT2,LSD.NET1,NET2, 10 

2NE T3 , NET 11 

INTEGER SEARCH, ACOMP,OKNOE 12 

EXTERNAL DIVE 13 

MCPF= 1 .0 1A 

MC PF 1= 1 #0 15 

ACOMP® 1 16 

LE VEL= 2 17 

R AE = 1 . 0A3 55 18 

RPE = 1 . 0A3 55 19 

RGE =6375, AA5 20 

GSE=.980665E-2 21 

PSI=0.0 22 

SE ARCH = 0 23 

KICK 1*1,0 2 A 

Kl CK2- 1.0 25 

KICK3=1.0 26 

TLT 1=0 ■ 0 27 

T LT2 = 0 .0 28 

TLT3=0 ,0 29 

I AE R0= 1 30 

NOE 1=1 31 

NOE 2 = 1 32 

NOE 3= 1 33 

M0DE=1 3A 

NE T 1 = -1 • 0 35 

NET2--1*0 36 

NET3=-1.0 37 

TN= 1,0 38 

PE N= • 1 39 

TLIMIT=1800.0 AO 

FELLOW=0.0 A1 

LEGT1=0.0 A2 

LE GT2 = 0 • 0 ^3 

LSD=0.0 AA 

6A5 TC=3600.0 A5 

GE= .980665E-2 A6 

VB = 0 *0 A 7 

DEGM.7A5329E-2 A8 

PY = 180 .0*DEG A 9 

TERR0R=1.0E-6 50 

IABC= 1 51 

OTBACK=2.0 52 

TIME I =0*0 53 

DELTI=A.O 5A 

K= 1.0 55 

WGM I N= 1 . 0F+ 10 56 

IOK 1= 2 57 

I0K2=2 58 

I0K3=2 59 

1 CALL INPUT (2HI1 ( I1,2H12,I2, 2H 13,13 ,2 HA 1 , A 1 , 2HA2 , A 2, 2HA3 , A3 , 2HB 1 ,B1 60 

1, AHNET1,NET1,AHNET2,NET2,AHNET3,NET3,3 61 

2HEW1,EW1, 3HEW2, E W2 , 3HE W3, E W3 , 3HET 1 ,ET 1 , 3HET2 , ET 2 , 3HET3 , E T 3 , 3HWEM , W 62 

3EM,3HWPL,WPL,5HV8C01,V8C01, 5HV 8C02 , V8C02 , 5HV8C03 ,V 8C03 , AHNOE 1 ,N0E1 63 

A,AHNOE2,NOE2, AHN0E3,N0E3,3HGSP,GSP,3HRGP,RGP,3HRPP,RPP,3HRAP,RAP,5 6A 

5HIAER0, I AERO, AHVTIT, VT I T , AHF W 1 2 , FW1 2 , AHMC PF , MC PF , 3HR AE , R AE , 3HRPE , R 6 5 

6PE , 3HRGE , RGE, 3HGSE , GSE , 5HT LTS I , T LT S 1 , 3HPS I , PS I , 5HDE LT I , DE LT I , 5HA ER 66 

70P , AEROP , 5HKICK1, KICK1, 5HK I CK2 , K I CK2 , 5HK ICK3,K ICK3 , AHT LT 1 , T LT 1 , AHT 67 

8LT2,TLT2, AHTLT3, TLT 3 , 6HSEA RCH , SEARCH , 5HMC PF l , MC PF 1 , 5HAC0MP , ACOMP , 6 68 

9HFELL0W, FELLOW, 3HL$D,LSD,5HLEGT1 ,LEGT1 ,5HLEGT2 ,LEGT2 , AHMODE , MODE , 3 69 

1HPEN,PEN, 6HTLIMIT.TLIMIT, 5HLEVEL , LEVEL ) 70 

GO TO ( 1112, 1111) , ACOMP 71 

1111 CALL INPUT! AHSIG1,SIG1, AHS IG2 , S I G2 ,AHS IG3 , S I G3 , 3HAT 1 , ATI, 3HAT2,AT2 72 

1 , 3HAT3, AT 3 , 3HB01 , B01 , 3HB02 » B02 , 3HB03 ,B03,3HCM1 ,CM1,3HCM2,CM2, 3HCM3 73 

2, CM3,5HTEMP1,TEMP1,5HTEMP2,TEMP2,5HTEMP3,TEMP3, 5HHE AT 1 , HE AT 1 , 5HHEA 7A 

3T2,HEAT2, 5HH EA T3 ,HEAT3,AHTEX1,TEX1 , AHTE X2 , TEX 2 , AHTEX 3 » T EX3 , 3HWK 1 , W 75 

AK1, 3HWK2, WK2, 3HWK3, WK3) 76 

1112 WRITE ( 6 , 100 )(VTIT(I),I=1,10) 77 

100 FORMAT! 1H1.20X, 10A6) 78 

IF (SEARCH. LE.O) GO TO 350 79 

IFINOEl.GT.il I 0K1= 1 80 

IF (N0E2.GT. 1 ) 1 OK 2= 1 81 

I F ( NOE 3.GT. 1 ) 1 0K3= 1 82 

350 WL=WPL 83 

MCFW1=(HL*MCPF1)-WL 8A 

WL=WL*MCPF 1 85 

WLSE1=LEGT1*FELL0W*LSD 86 

WLSE2=LEGT2*FELLOW#LSD 87 

WL=WL+WLSE1 88 

HP =W L 89 

NOE =N0E 1 90 

WE=EW1*FL0AT(N0E ) 91 

THRUST=ET1*FL0AT(N0E) 92 


ET=ET1 93 

E W= E Wl 94 

V8C0=V8C01 95 

A 12= A1 96 

112=11 97 

KICK=KICK1 98 

NET=NET1 99 

T LTSI = TLT1 100 

GS=GSP 101 

RG=RGP 102 

RP=RPP 103 

RA=RAP 104 

DELT=DELT I 105 

PSI2=PSI*0EG 106 

FH 12= THRUST/ (WL+WP+WE) 107 

DVGL=0.0 108 

GO TO (9999, 3333,4444,3333) .MODE 109 

3333 WG=WPL 110 

WEM=0.0 111 

GO TO 1113 112 

4444 I AER0=2 113 

AER0P=1.0 114 

WE M = 0 • 0 115 

WG=WPL 116 

GO TO 1113 117 

9999 IF(NET.LE.O.O)GO TO 25 118 

TN=FLOAT(NOE)/NET 119 

TN=TN+.9 120 

IOTN=TN 121 

TN= 10TN 122 

GO TO 9911 123 

25 TN= 1 .0 124 

9911 CALL 1NERT(WP,TANK,B0IS,SFMM) 125 

IF(NOE.EQ.l) GO TO 27 126 

WE =F LOA T ( NOE )*PEN*EW+ FLOAT ( NOE) *EW 127 

THRUST=FLOAT(NOE)*ET 128 

GO TO 21 129 

27 WE = E W 130 

THRUST=ET 131 

GO TO 21 132 

21 GO T0(221,66,22i) ,IABC 133 

221 WEWL=WE/WL 134 

TLTSR=TLTSI*DEG 135 

GM=GS*RG*RG 136 

VC=SQRT(GM/(RP*RG) ) 137 

AS=.5*IKICK*RP+RA) 138 

EP S= I R A-KI CK*RP ) / ( R A+K I CK*RP ) 139 

PS=AS*RG*(1.0-EPS*EPS) 140 

RLTS=PS/( 1 .0+EPS*COS I TLTSR) ) 141 

R0LTS=SQRT(GM/PS)*EPS*S1N( TLTSR) 142 

VSU=SQRT((GM/(AS*RG) >*(l.+2.*EPS*C0S (TLTSR)+EPS*EPS) / I l.-EP$*EPSU 143 

DVI=SQRT( (V8CO**2)+(2.0*GM/(RLTS )))-VSU 144 

87 CALL INERT(WP,TANK,BOI$,SFMM) 145 

DVIWGL=DVI+DVGL 146 

BMFST = EXP ( I -DVI WGL )/( I 12*GE ) ) 147 

BPROP = ( 1 • 0+ WE WL) /(BMFST*( 1 . 0+ A 12 )-A12) 148 

FW 12 = THRUS T/IBPROP*WL) 149 

WP2 = ( 1 • O-BMFS T ) • (BPROP*WL) 150 

IF ( ABS( 1.0-( WP/WP2) )-0.01) 88,88,89 151 

89 WP=WP2 152 

GO TO 87 153 

88 MFST * BMFST 154 

FW1 2= THRUST/ (BPROP*WL) 155 

2 TIME=TIMEI 156 

75 ALT0=FW12*GE 157 

XX3=GM/(PS**3) 158 

TH DT = SORT I XX 3 >* I I 1 .0+EPS*C0SI TLTSR ) ) **2 ) + ( VB/RLTS ) 159 

GO TO I 166, 566) , LEVEL 160 

166 CALL GLOSS 161 

TST=WP*I 12/THRUST 162 

DV IWGL=DVI +DVGL 163 

MFST=EXP(I-DVIWGL)/(I12*GE)> 164 

IF ( I ABC.EQ.2 ) MFS T= 1.0/MFST 165 

GO TO 60 166 

C FOLLOWING IS SETUP FOR INTEGRATION ROUTINE 167 

566 T ( 2 ) =TI ME 168 

T(3)=DELT 169 

T ( 4 ) =0*0 170 

T( 5 )=RLTS 171 

T ( 6 ) = THOT 172 

T ( 7 )=RDLTS 173 

T(8)=0.0 174 

CALL INTIT, 5,0, TERROR, 0,0,0 ,0,0, DIVE) 175 

55 CALL INTM 176 

533 VSQ=T(7)*T(7)+( (T< 5)*T(6) )**2> 177 

V=SQRT ( VSQ ) 178 

CAPE=.5*VSQ-GM/T( 5) 179 

H = T( 5 )*T{ 5 ) *T I 6) 180 

P=H»H/GM 181 

E=SQRTI 1.0+12. 0*CAPE*P/GM) ) 182 

GAMMA= ATAN (T(7)/(T(5)*T(6)) ) 183 

T0P=TI7)/(T(5)*T(6)) 184 
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BOT=1.0-(T( 5)/P) 

THETA=ARTN( TOP, BOT) 

IF (THETA.LT. O.O) THETA=THET A + 360 .0*DEG 
MF12=1.0-( ( FWl 2*T ( 2) ) / ( I 12 )) 

FWBO=FWI2/MF12 
IFJCAPEHO, 11,11 

10 V INF=0 . 0 
GO TO 50 

11 VIF2=2.0*CAPE 

V INF = SQRT ( V I F2 ) 

IF (VINF-V8C0)50,50,51 

50 H0LD2=T(2) 

H0LD4* T ( 4 ) 

H0LD5=T( 5) 

H0LD6=T ( 6 ) 

H0LD7 = T ( 7 ) 

HOLDB=T ( 8 ) 

GO TO 55 

51 T ( 2 ) =H0LD2 
T ( 3)=DTBACK 
T ( 4)=H0LD4 
T ( 5 ) =H0L05 
T ( 6 ) =H0LD6 
T ( 7 )=HOL07 
T ( 8 ) = H0LD8 
KS= 1 

CALL INT(T, 5, 1,0, 0,0, 0,0,0, DIVE) 

53 CALL I NTM 

VSQ=T (7 )*T( 7 )+( (T(5)+T(6))**2) 
V=SQRT(VSQ) 

CAPE=.5*VSQ-GM/TI 5) 

1 F ( CAPE. LT . 0 • 0) GO TO 53 
H=T( 5)*T( 5)*T(6) 

P=H*H/GM 

E=SQRT( 1.0+ ( 2. 0*C APE« P/GM ) ) 
GAMMA=ATAN(T(7)/(T(5)*T(6) )) 

TOP-T ( 7 ) /( T( 5) *Tt 6 ) ) 

BOT= 1 .0— ( T ( 5) /P) 

THETA=ARTN( TOP , BOT ) 

IF( THETA. LT. 0.0 ) THETA=TH6 TA+360 .0 *DEG 

REF=THETA-T( AJ-TLTSR 

A SS = 9Q ,0*DEG+AS I N ( 1.0/E) 

T0UT7=(ASS-REF)/DEG 

M F 1 2= 1 .0- ( (FW12*T( 2>)/( 112 )) 

VIF2=2.0*CAPE 

VINF=SQRT{ VIF2) 


185 

186 

187 

188 

189 

190 

191 

192 

193 

194 

195 

196 

197 

198 

199 

200 
201 
202 

203 

204 

205 

206 

207 

208 
2 09 
210 
211 
212 

213 

214 

215 

216 

217 

218 

219 

220 
221 
222 

223 

224 
2 25 
226 

227 

228 

229 

230 


IF(KS-5)54,54,56 
54 V IN8 ( KS ) = VI NF 
MF8(KS)=MF12 
TH0L8 ( K5 )=T( 2 ) 

T99 ( KS ) = T0l)T7 
KS=KS+ 1 
GO TO 53 

56 DO 57 I J- 1 » 4 

TVIN8( IJ)=VIN8(I J+l) 

TMF8( I J)=MF8< I J+l) 

TTH0L8 (IJ)=TH0L8(IJ+1) 

TT99( I J ) = T9 9 ( IJ+1) 

57 CONTINUE 

TV I N8 ( 5 ) = V I NF 
TMF8 ( 5 )=MF1 2 
TTH0L8 ( 5 )* T ( 2) 

TT99 ( 5 ) =T0UT7 
DO 58 14=1,5 
VIN8< 14)=TVIN8(I4) 

MF8 ( I4)=TMF8( 14) 

TH0L8U4)=TTH0L8U4) 

T 99 ( 14) = TT99 ( 14) 

58 CONTINUE 

DNN4=0 ,0 

IF(V1NF-V8C0)53,53,59 

59 CALL TAINT! VI N8 , MF 8 , V8C0, MFST , 5 , 2 , N I R , DNN4 , TH0L8 , TS T, T99 , TTTOT ) 

60 GO T0(61,62,61),IABC 

61 PROP=( 1.0+WEWL) /(MFST*( 1 . 0+A 1 2 ) - A1 2 ) 

WlGUS=PROP*WL 

FWTST=THRUST/HI GUS 
IF (ABS(FW12-FWTST)-.01)70,70,71 
71 FWl 2=FWTST 

WPH0= 1 .O-MFS T 

WO=PROP*WL 

WP=WPW0*W0 

CALL INERT ( WP , TANK, BO I $ , SFMM ) 

GO TO 75 

70 GO T0(777, 62,63) ,IABC 
777 WPWO= 1 • O-MFS T 
WTW0=WPW0*A12 
MUL= 1 .O/PROP 
W0=PR0P*WL 
WP=WPWO*WO 
WT=WTWO*W0 
WE=WEWL*WL 
VEM0S=V8C0/29.785 


231 
2 32 

233 

234 

235 

236 

237 

238 

239 

240 

241 

242 

243 

244 

245 

246 

247 

248 

249 

250 

251 

252 

253 
2 54 

255 

256 

257 

258 

259 

260 
261 
262 

263 

264 

265 

266 

267 

268 

269 

270 

271 
2 72 

273 

274 

275 

276 
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DV = I 12*GE*(-AI_0G< MFSTJ) 277 

TB= TST 278 

WG=WO 279 

VGL-OV-OVI 280 

I S AVE= 1 281 

IF(TB.LT.TLIMIT) GO TO 22 282 

NOE =NOE + 1 283 

I OK I s 2 284 

GO TO 9999 285 

22 GO TO! 110,666) , I0K1 286 

666 1 SAVE=2 287 

GO TO 176 288 

110 IF (WG-WGM1N) 176, 102, 102 289 

101 IF(NOE.LE.l) GO TO 111 290 

N0E=N0E-1 291 

WGM I N=WG 292 

GO TO 9999 293 

102 WGMIN=1.0E+10 294 

111 HR 1 TE ( 6 , 600 ) 295 

600 FORMATt 1H0,40X, 21H0EPARTURE FROM PLANET) 296 

WR I TE ( 6 , 60 1 ) 297 

601 FORMAT! 1H0, 3X , 7HENGINES, 3X,6HTHRUST,5X ,7HENG. WT ,4X ,4HV INF , 4X , 5HVG 298 

1L ,4X,5HDEL.V,4X,9HBURN T I ME ,4X , 3H I SP , 7X , 1HA , 7X , 3HMUL , 3X , 5HTANKS ) 2 99 

HRITE(6,602)OKNOE,OKTHR,OKHE,V8CO,OKVGL,OKOV,OKT8,I 12 , OKA 1 2 , OKMUl t 300 

10KTN 301 

602 FORMAT! 7X, I2,5X,F8. 1 , 3X, F8. 1 , 3X , F6 .3 » 2X , F6 .4 , 3X , F6.4, 3X »F8 . 2 , 5X , F6 302 

l.l f 3X,F5.3,3X,F6.4,3X,F6.1) 303 

WRITE(6,603) 304 

603 FORMAT! 1H0,3X,8HPAY L0AD,4X,6HW TANK,6X,6HW FU EL , 3X , 8HGR0SS WT,7X, 305 

11HP,7X,3HT/H,5X,8HMI0 FUEL , 5X ,7HL I FE WT , 3X , 6HTHETAP ) 306 

WRITE(6,604)OKWLtOKWT,OKWP,OKWGfOKP,OKFW,OKMFl ,OKHLStOKTT 307 

604 FORMAT! 3X , F9. 1, 3X, F8. 1,3X, F9. 1 , 3X , F9 . 1 , 3X , F6 . 3 , 3X , F6 . 4, 3X , F8 . 1 , 3X , 308 

1F7.1,3X,F8.1,3X,F6.2) 309 

IF! ACOMP.LE.l) GO TO 1113 310 

WRITE (6,605) OKTANK , OKBO , OK MM 311 

605 FORMATt 1H0 , 35H I NERT FRACTION A CONSISTS OF TANK= f F8.6 , 24H BOIL 0 312 

IFF + IN SUL AT ION=, F8. 6, 18H METERO ID SHI ELD=,F7 .5 ) 313 

1113 GO TO (700,800) , IAERO 314 

700 DELT=— 4.0 315 

WGMIN=1.0E+10 316 

PS I 2= 180 .0*DEG 317 

V 8C0 = V8C02 318 

N0E=N0E2 319 

ET=ET2 320 

EW=EH2 321 

WE=EW2*FL0ATf NOE) 322 

THRUST=ET2*FL0AT(N0E) 323 

NET=NET2 324 

A1 2= A2 325 

K I CK=K I CK2 326 

TLTSI=TLT2 327 

I 12-12 328 

I ABC= 2 329 

DTBACK=-2.0 330 

WL 1=WG+WEM 331 

HL=WG 332 

WP=2.0E5 333 

CALLINERT(HP,TANK,BOIS»SFMM) 334 

WTGUS=WP*A 12 335 

FW12=THRUST/( WL+WTGUS+WE+WEM) 336 

TLTSR=TLTSI*DEG 337 

GM=GS*RG*RG 338 

VC=SQRT(GM/(RP*RG) ) 339 

AS=.5*(KICK*RP+RA) 340 

EPS=(RA-KICK*RP)/!RA+KICK*RP) 341 

V P= SORT ! ( 2.0*GM/RG)*( ( 1. 0/RP) -( 1 .0/ (2 .0*AS ) ) 5) 342 

PS=AS*RG*( l.O-EPS+EPS) 343 

RDLTS=$QRT(GM/PS)*EPS*SIN( TLTSR) 344 

RLTS=PS/t 1.0+EPS*C0S(TLTSR) ) 345 

V SU= SORT ( (GM/(AS*RG) )*( 1 . +2. *E PS*COS (TLTSR ) +EPS*EPS ) / ( 1 .-E PS*EP S ) ) 346 

DVI=SQRT( (V8C0**2)+(2.0*GM/(RLTS )) )-VSU 347 

0VGL=0 .0 348 

143 CALL INERT(WP,TANK,BOIS,SFMM) 349 

WTGUS=HP*A12 350 

WE=EW2*FL0AT(N0E) 351 

THRUST=ET2* FLOAT (NOE) 352 

FW 12= THRUST/! WL+WTGUS+WE+WEM) 353 

DVIWGL=DVI+DVGL 354 

BMFST=EXP( (— OVIWGL ) /( I 12*GE) ) 355 

BPROP= ( 1.0+WE/WL1) /( BMFST* ( 1.0+A12 J-A12 ) 356 

WP=( 1.0-BMFST)*(BPROP*WL1) 357 

WTG 1 = WP» A12 358 

FW12=THRUST/!WL+WTG1+WE+WEM) 359 

IF! ABS ( 1.0-1 WTGUS/WTG1) )-. 01)141, 141,142 360 

142 WTGUS=WTG1 361 

GO TO 143 362 

141 WTGUS = A 1 2* WP 363 

M FS T= 1 . 0 /BMFST 364 

GO TO 9999 365 

66 WLC= WL+WTGUS+WE+WEM 366 

FW 1 2= THRUST/ WLC 367 

GO TO 2 368 


62 PROP=( 1.0+WE/WL1 J /((1.0/MFST)*! 1 . 0+ A1 2 ) -A 1 2 ) 369 

WBO=PROP*WL1 370 

WPC= ABS ( WBO-WLC ) 371 

WTC=WPC*A12 372 

IF (ABS I 1.0-WTC/WTGUS >-.01)65,65, 67 373 

67 CALL INERT! WPC, TANK, BOIS,SFMM) 374 

WTGUS=WTC 375 

GO TO 66 376 

65 PROP=( 1.0+WE/WL1) /( ( 1.0/MFST)*( 1 . 0+ A1 2 ) -A 1 2 > 377 

MUL = 1 .O/PROP 378 

WO=PROP»WL 1 379 

WP=WPC 380 

WT =WT C 38 1 

DV=1 12*GE*( ALOG(MFST) > 382 

VEMOS= V8C0/ 29 . 785 383 

WG=W0 384 

TB=ABS( TST) 385 

VGL= DV-DV I 386 

I SAVE = 3 387 

IF( TB.LT.TLIMIT) GO TO 23 388 

NOE=NOE+ 1 389 

I OK2=2 390 

GO TO 9999 391 

23 GO T0(210,667) , I OK 2 392 

667 I SAVE = 4 393 

GO TO 176 394 

210 IF(WG-WGMIN) 176,202,202 395 

201 WGM I N= WG 396 

IF (NOE .LE. 1 ) GO TO 211 397 

NOE =NOE- 1 398 

WTGUS=A12*( WL1-WEM) 399 

GO TO 9999 400 

202 WGMIN= 1.0E+10 401 

211 WRI TE (6, 508 ) 402 

WGMIN= 1.0E+10 403 

508 FORMAT( IH0,40X, 14HCAPTURE PHASE) 404 

WR I TE ( 6 » 60 1 ) 405 

WRITE (6,602)OKNOE,OKTHR,OKWE, V8C0, OKVGL ,OKDV ,OKTB , I 1 2 ,0K A1 2 , OKMUL , 406 

10KTN 407 

WRI TE (6, 613 ) 408 

613 FORMAT! 1H0,3X,8HPAY L0AD,4X,6HW TANK,6X,6HW FU EL , 3X , 8HGR0SS WT,7X, 409 

11HP»7X,3HT/W»5X, 6HE X MOD , 5X, 6HTHET AP ) 410 

WRI TE (6,604) OKWC, OK WT , OK WP, OKWG ,OKP ,OKFW ,WEM,OKTT 411 

I F( ACOMP.LE. 1 ) GO TO 987 412 

WRITE (6,605) OK TANK, OK80, OKMM 413 

987 I F ( MODE. EQ. 4) GO TO 645 414 

I ABC = 3 415 

WRITE(6,503) 416 

503 FORMAT< 1H0 , 39X, 20HDE PARTURE FROM EARTH) 417 

N0E=NOE3 418 

WE=EW3*FL0AT (NOE) 419 

THRUST=ET3*FLOAT(NOE) 420 

EW=EW3 421 

ET=ET3 422 

NET=NET3 423 

V 8C0= V8C03 424 

A1 2 = A 3 425 

KICK=KICK3 426 

TLTS I =TLT3 427 

112=13 428 

GS=GSE 429 

RG=RGE 430 

RP =RP E 431 

R A=RAE 432 

DELT=4.0 433 

P SI 2=0.0 434 

DTBACK= 2.0 435 

MCFW=(WG*MCPF)-WG 436 

WG=WG*MCPF 437 

W L=WG 438 

WL=W L+W LSE2 439 

W P = H L 440 

FW 12= THRUS T/IWL + WP+WE) 441 

DVGL=0 .0 442 

GO TO 9999 443 

800 I F ( MODE. EQ. 3) GO TO 987 444 

WR I T E ( 6 , 500 ) 445 

500 FORMAT! 1H0, 39X, 27HCAPTURE PHASE AERO BREAKING) 446 

PLW=WG+WEM 447 

PLW2=AER0P*PLW 448 

HSW=P L W2-P LW 449 

WRITE(6,501> 450 

501 FORMAT! 1H0, 3X,7H PAYLOAD, 5X,6HEX. MOD ,4X »8HHEAT SH . , 4X , 8HGR0SS WT ) 451 

WRITE(6,502 )PLW, WEM,HSW,PLW2 452 

502 FORMAT! 3X,F9. 1,3X,F8. 1, 3X,F8.1 ,3X,F9. 1 ) 453 

WG=P LW2 454 

GO TO 987 455 

63 WPW0=1.0-MFST 456 

WTWO=WPWO* A1 2 457 

MUL=1. O/PROP 458 

WO=PROP*WL 459 

WP=WPWO*WO 460 



WT=WTW0*W0 

461 


WE=WEWL*WL 

462 


V E MO S=V8 CO 729.785 

463 


DV= I 12*GE* ( -ALOG ( NFST ) ) 

464 


TB=TST 

465 


WG=*WO 

466 


VGL=*DV-DVI 

467 


I SA VE=5 

468 


IF(TB.LT.TLIMIT) GO TO 24 

469 


NOE =NOE+ 1 

470 


I0K3= 2 

471 


GO TO 9999 

472 

24 

GO TO (310,668),I0K3 

473 

668 

I SAVE=6 

474 


GO TO 176 

475 

310 

IF (WG-WGMINI176, 302,302 

476 

301 

IF (NOE.LE.l) GO TO 311 

477 


N0E=N0E~1 

478 


WGM IN=WG 

479 


GO TO 9999 

480 

302 

NOE =NOE+ 1 

481 


I OK 3= 2 

482 


‘HGM IN = 1 .0E + 10 

483 


GO TO 9999 

484 

311 

WR 1 TE < 6 , 60 1 ) 

485 


WGM IN= 1 ,0E+ 10 

486 


WRITE (6,602>OKNOE,OKTHR,OKWE,V8CO,OKVGL,OKDV,OKTB, 1 1 2 ,0KA1 2 ,OKMUL , 

487 


10KTN 

488 


WR I TE ( 6 , 62 3) 

489 

623 

FORMATt 1H0, 3X,8HPAY L0AD,4X,6HW TANK,6X,6HW FUEL , 3X ,8HGR0SS WT,7X, 

490 


11HP,7X,3HT/W,5X,8HMI0 FUEL, 5X , 7HL I F E WT , 5X .6HTHETAP ) 

491 


WRITE(6*604)0KWL» OKWT ,OKWP,OKWG ,OKP ,0KFW , 0KHF2 * OK WL2, OKTT 

492 


IF (ACOMP.LE. 1) GO TO 99 

493 


WRITE (6,605 I OK TANK , OKBO, OKHM 

494 

99 

GO TO 645 

495 

176 

OKNOE=NOE 

496 


OKTHR=THRUST 

497 


OKWE-WE 

498 


DKVGL=VGL 

499 


OKDV-OV 

500 


OKTB=TB 

501 


OKA 1 2= A12 

502 


OKMUL-MUL 

503 


OKTN=TN 

504 


OKWL=WL 

0KWC=WL1 

505 


OKWT=WT 

507 


OKWP=WP 

508 


OKWG=WG 

509 


OKP=PROP 

510 


OKFW=THRUST/WG 

511 


OKMF 1=MCF W 1 

512 


0KWLS=WLSE1 

513 


OKTT=TTTOT 

514 


OKTANK=TANK 

515 


OKBO=BOIS 

516 


OKMM=SFMM 

517 


0KMF2=MCFW 

518 


0KWL2= WLSE2 

519 


GO TO ( 101, 111,201,211,301,311), ISAVE 

520 


END 

521 


C APPROXIMATE GRAVITY LOSS ROUTINE 

SUBROUTINE GLOSS 
C 

C THIS ROUTINE IS AN APPROXIMATION FOR THE GRAVITY LOSSES, ANO IS 

C USED ONLY FOR L EV EL= 1 . WHEN L E VEL =2 THE GRAVITY LOSSES ARE SET TO 

C 0 AND RETURNED TO THE MAIN PROGRAM WHERE THE TRUE GRAVITY LOSSES 

C ARE FOUND BY INTEGRATION OF THE EQUATIONS OF MOTION. THIS ROUTINE 

C MAY BE REPLACED BY THE USER WITH HIS OWN APPROXIMATION IF DESIRED, 

C TO DO THIS ALL NEEDED VALUES MUST BE PLACED IN THE COMMON BLOCK 

C / BO XGL / HERE AND IN THE MAIN PROGRAM THE VALUE OF THE GRAVITY 

C LOSSES MUST BE PLACED IN DVGL BEFORE RETURN TO THE MAIN ROUTINE 

C 

COMMON ALTO, P SI I, I 1 2 , GE , GM , A LT , T ( 70 > , PS I 2 , K ,GS 

COMMON /BO XGL/ DVI , FW12, V8C0,VSU ,DVGL , L EV EL , T ST »RG,RLTS»AS 

REAL 312 

GO TO (2,5), LEVEL 

2 DVGl=.036+.0693*( ( V8C0/ VSU) **2 ) 

FACTl=1.0/( AS*AS) 

FACT2= (GS/GE )**2 
P 1 = FW12*.5 
P2 = I 12*GE 

P5=0VG1*DVG1*VSU*FACT1*FACT2/(P1*P1) 

P6 =E XP { DV I /P 2 ) 

G=DVGL 

9 P9=1.0 + (P6*EXP(G/P2) ) 

PART 1 = 2. 0*P5*P6*EXP< G/P2)/P2 
PART2=P9**3 
F0FG=G-(P5/(P9*P9) ) 

FPRI M= 1.0+ (PARTI/ PART2 ) 

G1=G-(F0FG/FPRIM) 

IF ( ABS( 1.0-G1/GI-.001 ) 25,25,26 
26 G=G 1 

GO TO 9 

25 DVGL= . 5* ( G 1+G ) 

5 RETURN 

END 
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C THIS IS THE DERIVATIVE ROUTINE FOR THE INTEGRATION ROUTINE 

SUBROUTINE DIVE 

COMMON ALTO , P S I 1 , I 1 2 , GE , GM , A L T , T ( 70 ) , PS 1 2 »K , GS 
DEG=. 01745329 
REAL K, I 12 
TOP = T< 5)*T<6) 

BOT= T ( 7 ) 

SI T= AR TN ( TOP » BOT ) 

IF(SIT.LT.O.O) SIT=SIT+360.0*DEG 
PSI 1=PSI2+K*SIT 
PRT 1 = T ( 2 ) / ( I 12*GE) 

PRT2=1. 0/ALTO 
ALT=1.0/(PRT2-PRT1) 

XX 1= T ( 5)*T(6) + T(6)-< GM/(T( 5 J*T( 5) ) > 

XX2=ALT*C0S(PSI 1 ) 

T ( 12)=XX1+XX2 

T( 11) = (ALT*SIN( PSI1)/T(5))-(2.0*T(7I*T(6) )/T(5) 

T t 9 ) = T ( 6 ) 

T( 10 > = T I 7 ) 

T ( 131=0.0 

RETURN 

END 
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C INERT SUBROUTINE 38 

SUBROUTINE I NERTI WPI » TANK,BOIS ,SFMM ) 3q 

C 40 

C SUBROUTINE INERT COMPUTES THE INERT FRACTION FOR TANKAGE, BOIL-OFF 41 

C AND INSULATION, AND METEROID PROTECTION (INCLUDES WHIPPLE BUMPER) 42 

C EFFECT. AC0MP=1 INERT FRACTION IS INPUT CONSTANT AC0MP=2 INERT 43 

C FRACTION IS COMPUTED. IABC IS CONTROL FOR PHASE OF MISSION BEING 44 

, C COMPUTED 1=PLANET D6 PAR TURE , 2 = PL ANE T CAPTURE , 3=EARTH DEPARTURE 45 

C 46 

! COMMON /B0X1/A1,A2»A3,A12» IABC , ACOMP , S IG1 ,S IG2 , SIG3, AT I , B01 ,CM1 , AT2 4 7 

1, B02,CM2,AT3,B03,CM3,TEMPl f TEMP2,TEMP3,HEATl,HEAT2,H6AT3,TEXl,TEX2 48 

I 2»TEX3»WKl» WK2, WK3 » TN 49 

INTEGER ACOMP 50 

WPI=WPI*. 4535924 51 

GO TO (1,21, ACOMP 52 

1 GO TO (3, 4,5), IABC 53 

3 A 12=A1 54 

GO TO 90 55 

4 A 12=A2 56 

GO TO 90 57 

5 A 12= A3 58 

GO TO 90 59 

2 GO TO (7, 8, 9), I ABC 60 

7 51 G=S IG 1 61 

i A=AT1 62 

| B= BO 1 63 

1 C=CM1 6 4 

| TEMP=TEMP1 65 

EL=HEAT1 6fe 

! £XTM£=TEX1 67 

| WK=WK1 68 

1 GO TO 80 69 

8 S 1G=SIG2 70 

I A=AT2 71 

B=B02 72 

C=CM2 73 

TEMP=TEMP2 74 

EL=HEAT2 75 

EXTME=TEX2 76 

WK=WK2 77 

GO TO 80 78 

9 S I G=S IG3 79 

A=AT3 80 

B= BO 3 81 

C=CM3 B2 

TEMP=TEMP3 83 

EL=HEAT3 84 

EXTME=TEX3 85 

WK=WK3 86 

80 TAU=A/( SIG**.533) 87 

83 TANK=< <TAU) *< TN**, 1) )/<WPI**. 1 ) + < WK*TN ) /WP I 88 

81 B0T=(WPI**.33333)*( S1G**. 666667) 89 

BOIS=(B/BOT)*SQRT(EXTME*TEMP/EL) 90 

T0P=(EXTME**.25)/( SIG**( 5.0/6. 0) ) 91 

T0Pl=1.0/( WPI** ( 1.0/6. 0) ) 92 

SFMM=C*T0P*T0P1 93 

I F ( SF MM- TANK ) 30 , 30 , 31 94 

30 SFMM=0.0 95 

GO TO 91 96 

31 WTEST=1.333333*TANK 97 

IF( SFMM-WTEST) 32,32,33 98 

32 SFMM® SF MM- TANK 99 

GO TO 91 100 

33 WTEST2=( 16.0/3.0) *TANK 101 

IF( SFMM-WTEST2) 34, 34, 35 102 

34 SFMM=TANK/3. 0 103 

GO TO 91 104 

35 SFMM=( SFMM/4.0)-TANK 105 

91 A12=T ANK+BOI S+SFMM 106 

WWWP=WPI/. 4535924 107 

90 WPI=WPI/. 4535924 108 

110 RETURN 109 

END HO 
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